Tight-binding models for the new iron based superconductor materials 
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The rich novel materials class of iron based superconductors turned out to exhibit a very complex 
electronic structure, despite of the simplicity of their crystal structures. For various approaches 
to study the instability against magnetic order or superconductivity, a real space description of 
the electronic structure is required. Here, the bonding situation and the orbital structure of the 
electronic state are analyzed and minimum tight-binding models quantitatively correctly describing 
the low-energy electronic structure are provided. 
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I. INTRODUCTION 

The discovery of high-temperature superconductivity 
in La(0,F)FcAs — immediately raised great expectations 
to advance the theory of high-temperature superconduc- 
tivity in general, since the crystal structure appeared 
quite simple and on a first glance displayed similarity 
to that of the cuprates. It quickly turned out, however, 
that there are also striking differences, in particular of 
the correlated nature of the electron state. The analysis 
of many problems in this respect demands real space ap- 
proaches based on tight-binding (tb) models. An early 
paper in this direction used a two-orbital tb model.— The 
need for gradually improving such models was already 
clearly seen in this paper. Unfortunately, improvements 
turned out to grow terribly in complexity, and hopping 
matrix elements up to the fifth neighbor matter^ The sit- 
uation clearly calls for a detailed analysis of the electronic 
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structure in the tb language as complete as possible. The 
present paper tries to contribute to this task. 

Among the new iron based superconducting materials 
there are essentially four structural families christened 
11, 111, 11 11 j and 122 according to their stoichiomctry 
in the undoped cases. FeSe and FcTe are representa- 
tives of the first family,— LiFeAs and NaFeAs represent 
the second^ REOFePN, where RE is Y, La or a light 
rare earth element and PN is P, As or Sb, belong to the 
third,— & and AFe2As2, where A is Ca, Sr or Ba, belong 
to the last.— All families exhibit a tetragonal structure 
at room temperature which slightly distorts below the 
Neel temperature of the anti-ferromagnetic order. The 
first three families have the non-symmorphic space group 
P4/nmm in common, while the last family crystallizes in 
the symmorphic space group I4/mmm. The unit cells of 
the first three families as well as the tetragonal cell of the 
last (which contains two unit cells) are shown in Fig. [TJ 





Figure 1: From left to right: unit cells of FeSe, LiFeAs, LaOFeAs and BaFe2AS2 (the tetragonal cell contains two 

unit cells in the latter case). 



The generic structural element is a metallic iron chalco- plane sandwiched between two chalcogen/pnictogen 
genidc/pnictide layer consisting of a square iron atom (c/p) atom planes in such a way that every Fe atom 
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is in the center of a c/p tetrahedron. The families 11 
and 111 consist of neutral metallic layers only, stacked 
on top of each other in the direction of the tetragonal 
axis and only weakly bonded in this direction. Thereby 
LiAs is isoelectronic with Se. They should easily cleave in 
planes perpendicular to the tetragonal axis. In contrast, 
in the 1111 and 122 families the FcAs layer is charged 
and anionic and is intercalated with a non-conducting 
cationic layer, e.g. LaO in the 1111 family and Ba in 
the 122 family. Hence, these materials cannot neutrally 
cleave between intact atom layers. 

Phenomenologically, superconductivity seems to be 
quite robust in the (doped) families 1111 and 122, but 
more sensitive to imperfections in the families 11 and 
111. This at least goes in line with a remarkable struc- 
tural difference revealed in Table |TJ For representatives 
of the four families this table shows the nearest in-plane 
Fe-Fe distance d in A, and the ratio r of the distance of 
c/p layers from the Fe layer related to the in-plane Fe- 
Fe distance. This latter ratio would be 0.5 in the case 
of regular c/p tetrahedra. As is evident from Table [H 
the cationic interlayers of the 1111 and 122 families put 
the Fc layer under considerable tensile strain, and in line 
with that the c/p tetrahedra are considerably contracted 
in the direction of the tetragonal axis. The ratio r is by 
the way unusually badly reproduced by lattice relaxation 
within non-magnetic density functional theory (DFT), 
and it seems to be somewhat sensitive to magnetic or- 
der. Which of the Fe 3d bands cross the Fermi level and 
the sequence in energy of these bands depends on the 
actual value of r and often differs for the experimental 
value and for the DFT relaxed value. 



family 


representative 


d 


r 


11 


FeSe 


2.6653 


0.5350 


111 


LiFeAs 


2.6809 


0.5614 


1111 


LaOFeAs 


2.8497 


0.4620 


122 


BaFe 2 As 2 


2.8019 


0.4855 



Table I: Structural parameters of the four families as 
described in the text. 

All results presented in this paper refer to the experi- 
mental structure data and to the non-magnetic state. 

The aim is to provide tb models of the low-energy back- 
ground electronic structure on the basis of which collec- 
tive behavior like magnetism and superconductivity may 
develop and may be analyzed by many-body theories. It 
will be shown that the low-energy physics is essentially 
more two-dimensional (2D) in the families 1111 and 122 
than in the other two, again distinguishing the lower two 
families of Table [J from the upper two. 

On Fig. [5] the essential atomic-shell resolved partial 
densities of states (DOS) together with the total DOS are 
shown for the experimental structure parameters of FeSe 
(a = 3.7734 A, c = 5.5258 A, z Se = 0.267, z Fe = 0)4 All 
band structures in this work are obtained with the high 
precision full potential local orbital code FPL082 using 



the Perdew-Wang 92iH version of the local density ap- 
proximation (LDA) density functional. Generalized gra- 
dient approximation (GGA) results would somewhat dif- 
fer in the lattice relaxation by total energy minimization. 
Since, however, for given fixed structural parameters the 
results considered here hardly differ at all, and all results 
presented in this text are for the experimental structures, 
LDA was used throughout. 

For what follows it is relevant that the all electron ap- 
proach FPL08 uses a carefully optimized 'chemical' basis 
(one basis orbital per atomic core or valence orbital) plus 
at most one polarization orbital per band. The quality 
of the chemical basis results in occupation numbers of 
polarization orbitals usually well below 0.02 (per spin), 
in most cases even much lower (without loss in accuracy 
which competes with any of the advanced all-electron 
full-potential tools presently available. For instance, in a 
recently published workii results for the 1111 compound 
LaOFeAs obtained with different band structure meth- 
ods were compared and very good aggreement between 
FPLO and FPLAPW results was assessed). This is the 
reason why projections on the FPLO chemical basis or- 
bitals have a lot of chemical relevance, and also produce 
nearly optimally localized Wannier orbitals from band 
structures in an approach described below in Sec. IV. 

The situation shown in Fig. [5] which is generic for all 
four families considered here can be understood as fol- 
lows: Since there are two Fe atoms per unit cell, there 
is Fe-Fe homo-covalency. Put the structural layers in 
the icy-plane stacked in z-direction. For zero wave vector 
there are in-plane bonding combinations of Fe-3c? orbitals 
with angular dependence xy/r 2 and (x 2 — y 2 )/r 2 , respec- 
tively, with the same coefficient on both sites and anti- 
bonding combinations with sign-alternating coefficients. 
For the angular dependence xz/r 2 and yz/r 2 it is just the 
other way round, while the orbitals ~ (z 2 /r 2 — 3) form es- 
sentially non-bonding bands. Hence, there are four Fe-Fe 
bonding bands, four Fe-Fe anti-bonding bands and two 
non-bonding bands. The Se-4p orbitals may couple with 
the Fe-Fe anti-bonding combinations only, developing Fe- 
Se hetero-covalency with predominantly Se-4p bands be- 
low —3 eV and predominantly Fe-3d bands above the 
Fermi level (the latter put at zero band energy). The 
Fermi level itself falls in a Fe-Fe covalency pseudo-gap 
between the bonding plus non-bonding bands of essen- 
tially pure Fe-3d character and the anti-bonding bands 
of hybridized Fe-3d-Se-4p character. This picture pre- 
vails in large parts of the Brillouin zone (BZ) and only 
modifies in the vicinity of the BZ-edges parallel to the 
z-axis (through point M of Fig. [5] below). It results in 
rather small hole Fermi surfaces (FS) of Fe-Fe bonding 
bands around the fc z -axis and in likewise small electron 
FS around the BZ edges parallel to this axis, and hence 
in the pseudo-gap at the Fermi level. The materials are 
not far from semi-metals. 

There is also Fe-4s-Se-4p covalency with a gap of about 
7 cV as seen on the bottom of Fig. [5] The modifica- 
tions in the other families compared to 11 are gradual, 
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Figure 2: Partial and total DOS for FeSe. The Fermi 
level is at eV. 
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Figure 4: Partial and total DOS for LaOFeAs. 
Experimental structure parameters from RefJ^. 



see Figs. [3H51 Since the As-4p states are about 1 eV 
higher than the Se-4p states, the gap below —2 eV is 
less pronounced or absent in the arsenides. In LiFeAs 
the pseudo-gap at the Fermi level is less pronounced too 
as the Fermi radii are somewhat larger (larger energy 
overlap between Fe-Fe bonding and anti-bonding band 
groups). In LaOFeAs and BaFe2As 2 there is also no gap 
above 2 eV as the La-4/ states and the Ba-5d states, re- 
spectively, come down close to the Fermi level. (About 
0.4 Ba-5d states are occupied in BaFe2As2 while, when 
doping with K, the K-3d states remain empty; this is why 
one K atom when replacing Ba donates only about 0.6 
holes, or 0.3 holes per Fe atom, into the Fe-3d bands.) 
Moreover, the 0-2p states of LaOFeAs appear below 
—2 eV in addition to the As-4p states. In the range 
between —2 eV and the Fermi level there is not so much 
difference in the bands being almost purely Fe-3d in all 
cases. 
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energy in eV 

Figure 3: Partial and total DOS for LiFeAs. 
Experimental structure parameters from Refj^. 

The whole situation with Fe-3d-Fe-3<i and Fe-3<i-c/p- 
4s covalency is illustrated for the generic example of FeSe 
in Fig. [SI The projection of the LDA bands on the FPLO 
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Figure 5: Partial and total DOS for BaFe2As2. 
Experimental structure parameters from Rcf. 14 . 



chemical basis orbitals is given as line thickness in color 
code. (The diameter of the color balls in the legend of 
the figure indicates the thickness for 100 percent orbital 
content in the band wave function.) The cyan (c/p-pz) 
band on the line T-Z which disperses down close to the 
Fermi level has a black component (Fe-xy) hidden be- 
hind the cyan and amounting to about half of the band 
state near Z. The narrowness below the Fermi level of the 
non-bonding Fe-z 2 bands is clearly seen. Again seen is 
that the Fe-Fe bonding bands below the Fermi level do 
not hybridize with the c/p-4p orbitals while Fe-Fe anti- 
bonding bands strongly hybridize. Due to the far extend 
of the radial orbitals, the c/p-4p states extend even above 
the Fe-3d bands, over an energy range of more than 10 
eV. It is also seen that the Fermi level is crossed by Fe 
xy, xz and yz bands only with a small admixture of c/p 
Apz orbitals. Right below Fermi level Fe x 2 — y 2 and z 2 
admixture sets in, so that one fails to reproduce Fermi 
velocities even qualitatively without consideration of the 
latter orbitals and their corresponding bands. 

Remarkable in this respect is also the orbital 
anisotropy of the FS, again in stark contrast to the case 
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Figure 6: Orbital projection of the LDA band structure 
of FeSe. (See Fig. M for the points T, X, M; Z, R, A are 
above the former at k z = ir.) 
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Figure 7: Partial and total DOS for a hypothetical 
structure obtained by removing the Se atoms from FeSe. 



of the cuprates. While the three FS hole cylinders around 
T-Z appear more or less orbitally isotropic with one of the 
outer cylinders and the inner cylinder of xz — yz charac- 
ter (of course with their mixing rotating together with k) 
and the other outer cylinder of xy character, the outer 
FS electron cylinder around M in the k z = plane ap- 
pears totally orbitally anisotropic with xz — yz character 
in the directions towards X and nearly pure xy character 
in the directions towards T. The inner electron cylinder 
is again more or less orbitally isotropic with xz — yz char- 
acter. In the k z = tt plane the inner and outer electron 
cylinders have interchanged their orbital character. 

In all those considerations the presence of the c/p 
atoms is essential. Without there presence, the unit cell 
would reduce to containing one Fe atom only, and no 
interband Fe-Fe covalency could develop. In Fig. [7] the 
DOS of a hypothetical lattice is shown with the Se atoms 
removed from the FeSe crystal, but the Fe atoms left 
at their positions. The Fc-3g? occupation remains nearly 
the same (Mullikan analysis with the FPL08 basis re- 
sults in charge transfer between Fe and c/p of less than 
0.1 electron charge/atom), but no pseudo-gap whatso- 
ever appears as there are now only five Fe-3d bands in 
the doubled BZ. There is also no covalency band gap of 
the Fe-4s states any more. Fe-Fe covalency appears only 
due to the doubling of the unit cell caused by the pres- 
ence of the c/p atoms, due to their additional potential 
and due to hybridization with c/p-p orbitals. Also, the 
width and depth of the pseudo-gap are quite sensitive 
to the c/p position, that is, to the actual value of the 
Wyckoff parameter z c / p . 

In a vicinity of half an eV around the Fermi level, the 
band structure in the 1111 family exhibits a pronounced 
quasi-2D character. It is somewhat less 2D and in a nar- 
rower energy window at the Fermi level in the 122 family 
and much more 3D in the 11 and 111 families. 

In the next section the implications of symmetry on 
the electronic structure arc analyzed in detail. This pro- 
vides some key to the construction of reduced tb models. 
In Sec. Ill, a minimum 2D tb model is derived for reason- 



able quantitative approximations of the low-energy band 
structures of the Fe-c/p layers (the focus being essentially 
on excitation energies below 0.1 eV), and parameters are 
given for representatives of all four families. These are 
based on a Wannier function representation introduced 
in Sec. IV. The families 11 and 111 exhibit quite sizable 
3D dispersion close to the Fermi level. Corresponding 3D 
tb models for these cases are given in Sec. V, and finally 
the results are summarized in Sec. VI. 



II. SYMMETRY OF THE IRON 
CHALCOGENIDE/PNICTIDE LAYER 

The 11 family consists of the Fe-c/p layers only, generic 
for all four families. To be specific, FeSe is taken as 
example in what follows. We start with a more detailed 
description of the FeSe structure: 

The two Fe and Se positions per unit cell are symboli- 
cally distinguished by a sign. Atom positions are (Fig. [5]) 

S±Fe = 7 (=FHi±-Ra), 

i w 

S±Se = j(=F-Rl T R2) ± ZSeR-3- 

The vectors T x and T y of Fig. [5] are not lattice vectors 
of the full structure. However, in combination with a re- 
flection n on the (a;, y)-plane they nevertheless form non- 
symmorphic symmetry elements (n\T x ), (n\T y ) (diagonal 
glide plane, (n\T x )(j)(r) = (p(n~ 1 r — T x ), n^ 1 = n). Ap- 
plied to a Bloch state ipkv{f), any space group element 
yields (A\T A +Ri)^ k „(r) = exp(i(Ak)-Ri)^' Ak (r), where 
v is the band index, and ipAk( r ) ' 1S a linear combina- 
tion of states ipAkv' (f) with the restriction EAkv' = £kv 
on the band energies. (T„ = T x , some one-one rela- 
tion A <-> Ta must be chosen, (n\T y ) = (n\T n + R 2 ). 
Moreover, the result of (j4|Ta + Ri)ipku( r ) is based on 
a chosen phase relation between ipk^if) and ipAk( r )' 
which choice is made here for the sake of simplicity 
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Figure 8: Tetragonal plane of nine unit cells (full-line 
quares) spanned by lattice basis vectors i?i and R2 of 
FeSe. The open circles mark Se positions above and the 
shadowed circles below the drawing plane. The Fe sites 
(full circles) are located in the drawing plane. If the Se 
sites are neglected, the half-size unit cell is spanned by 
T x and T y as shown in the left lower corner. 



different from that later in (jlQ|) and (fT6]k the results 
considered here do not depend on this choice.) One 
has (A\T A )(B\T B ) = (AB\T AB + Rab), where R AB 
is a lattice vector. If, for the wave vector group of k, 
{(A\Ta + Ri) I Ak — k + Ga} where Ga is a reciprocal 
lattice vector, one defines the multiplier representation 
pt, v {k) = exp(-ifc • Ri){rl) hvl \{A\T A + Ri)\ip h v), then 



P A P B 



exp(ifc ■ Rab) P 



AB 



(2) 



follows. Time inversion T is defined as the anti-unitary 
operator T<p{r) = <jf(r) and hence Tipk v = ip-kTv and 
T(A\T A + Ri) = (A\T A + Ri)T. Now, if m 2 is the re- 
flection plane through iron sites sending R2 — ► —R2, 
then the anti-unitary operator = Tm2{n\T n ) is an 
element of the wave vector group for any wave vector 
k with k ■ Ri = ±7r, that is, on the face of the Bril- 
louin zone (BZ) with central point X (Fig. [5] below), and 
p & (k)p & (k) = —1 for any such k since O 2 = (1\R\). 
The existence of an anti-unitary operator with square 
equal to —1 which commutes with the Hamiltonian is 
an analogue to Kramers degeneracy^ and an immediate 
consequence here is a twofold degeneracy of all bands on 
the four side faces of the BZ. (Ri and R 2 may be inter- 
changed in the above considerations, with m 2 replaced 
with nil.) One also has — {mi{n\T n ) + (n\T n )nii)/2 = 
2\, (m 2 (n\T n ) + {n\T n )m 2 )/2 = 2\ where 2\ is a screw 
rotation by 7r around the axis through the lattice vector 



R2 of Fig. [U accompanied by a glide by the vector R2/2, 
and 2\ is the corresponding screw rotation around R\ as 
given in that figure. The notation of space group ele- 
ments follows Hermann-Mauguini£. We prefer to use the 
two reflections in the considerations of hopping matrix 
elements below. 

This is all generally true for the P4/nmm space group 
(which is also the space group of a single FeSe layer of 
Fig. [8]). Hence, it holds for the 11, 111 and 1111 materials 
(FeSe, LiFeAs and LaOFeAs as examples). It does not 
hold for the 122 materials (e.g. BaFe2As2) which crys- 
tallize in the I4/mmm structure (body-centered tetrag- 
onal). The space group I4/mmm is symmorphic, and 
hence there appear no multipliers different from unity in 
representations of wave vector groups and no two-fold 
degeneracies on whole BZ faces. These degeneracies only 
reappear, if one neglects z-axis dispersion. In the follow- 
ing we first focus on the layer structure of Fig. [5J which 
is representative for all four families, if z-axis dispersion 
is neglected. 

The full crystal potential V must be invariant under all 
symmetry operations including lattice translations and 
glide reflections. With respect to the latter, it can be 
decomposed into a part V+ that is invariant under both 
the T x - or T^-translations and the reflection n on the 
(x, y)-plane separately and a part VL that is alternating 
with respect to both separate operations. This potential 
representation will be used in the next section. 

Further important symmetry transformations of all 
structures with the P4/nmm space group are a 7r-rotation 
2 X around the a;-axis of Fig. [51 a 7r-rotation 2 V around a 
parallel to the y-axis through Fe sites and the spatial 
inversion 1 with its center halfway between any of the 
pairs ±X of atoms (all pairs have the same midpoint). 
Non-zero Hamiltonian matrix elements must either be in- 
variant under symmetry transformations of the structure 
or transform into each other. 

Due to the presence of the glide plane (n\T n ), all atom 
positions within the unit cell of the layer of Fig.[5]or more 
generally of a crystal with P4/nmm space group come 
in symmetry equivalent pairs ±A like ±Fe and ±Se in 
Fig. [51 Denote local site orbitals of a local orbital basis 
correspondingly by a sign superscript as ± . (This super- 
script is to be well distinguished from the sign subscript 
on V± for the component of the crystal potential which 
has a different meaning defined above, although related 
to the same symmetry element.) Expanded in this ba- 
sis, the Hamiltonian may be brought into a block matrix 
structure 



H 




(3) 



if one denotes hopping matrix elements by 



where 

for hopping a distance rT a 



sT., 



ii?3, r, s integer, be- 
tween local site orbitals <f>i and <f>j (of either superscript), 
then hoppings with r + s even (accompanied with corre- 
sponding Bloch phase factors) enter H ++ and H~~ while 
hoppings with r + s odd enter and H '". 
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Note that the anti-unitary transformation Tl is an el- 
ement of the wave vector group for every wave vector fe, 
this time with (p T1 ) 2 = +1 and hence not causing a de- 
generacy. The local site orbitals may always be chosen 
so that 



There are 4 distinct onsitc Hamiltonian diagonal ma- 
trix elements (the subscript i being again the number 
from ©) 



£4 = £3, 



ti+5 



(8) 



Tl</. d 



(4) 



this time also for the 122 structure, in all cases with re- 
spect to the center of inversion 1 half-way between near- 
est Fe neighbors. 

For instance, if in a natural way the local site or- 
bitals are parity eigenstates, appropriate imaginary fac- 
tors at odd parity orbitals do. Then, for the hopping ma- 
trix elements one has {(jr\H\(jr) = (Tl(j>-\H\Tl<f>-)* = 
{(j}+\H\<j) + )* , and {<jr\H\<i>+) = {<p + \H\ </>")*. It immedi- 
ately follows that with this basis choice (which is always 
possible and even the natural one) 



H 




(5) 



Since generally H + = H + , it also follows that the 
block matrix must be symmetric: 



H 



H- 



(6) 



where H l means the transposed of H . This reduces the 
calculation of Hamiltonian (and overlap) matrix elements 
by a factor of two. 



III. TIGHT-BINDING BANDS FOR THE IRON 
CHALCOGENIDE/PNICTIDE LAYER 

In this section tb models for the iron c/p layers are 
considered. Again, to be specific, the discussion is first 
for FeSe and then generalized. 

Fe 3d-orbitals arc denoted by their angular character 
and by an upper sign indicating centering at site S±F e 
as in the previous section for the general case. For in- 
stance (ixz) + means a real orbital with angular depen- 
dence xz/r 2 centered at S+Fe, multiplied with the imag- 
inary unit factor later being used to get real Hamiltonian 
matrices; (z 2 ) abbreviates the real orbital with angular 
dependence (3z 2 — r 2 )/r 2 . The following site orbitals </> ± 
within a unit cell are introduced which obey (|4"|): 



1 
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(xy) , 


2 


(x 2 -y 2 )+, 


7 


(x 2 - v 2 r 


3 


(ixz) + , 


8 


(—ixz)~, 


4 


(iyz)+, 


9 


(-iyz)~. 


5 


(^ 2 ) + , 


10 


(z 2 )-. 



(7) 



The numbering of orbitals will later be used for (lower) 
matrix indices of the hopping matrices tij. If in what 
follows no superscript is applied, then means an orbital 
either on site S+Fe or on S-Fe- 



while due to orbital onsidc orthogonality there are no 
onside Hamiltonian off-diagonal matrix elements. The 
10 x 10 Hamiltonian matrix in site orbital representation 
has an obvious 5x5 block structure [Gj . 

Since only one layer is considered, only in-plane hop- 
pings t^j a distance rT x + sT y figure, a fact relevant in 
what follows. Now one has to distinguish matrix ele- 
ments, the orbital products of which arc even with re- 
spect to reflection on the (x, y)-plane and which hence 
are coupled by the kinetic energy operator and by V+, 
from the ones, the orbital products of which are odd 
with respect to reflection on the (x, y)-plane and which 
hence are coupled by V—. The former are (superscripts 
dropped, i,j = 1, . . . , 5) tu, h 2 , *i5, *25, *34, their trans- 
posed and the corresponding ti+sj+s, U,j+5 and ij+sj. 
The later are Uj, i = 1,2,5; j = 3,4 their transposed 
and the corresponding ^+5^+5, tij+5 and ij+5,j. 

For the former one has immediately ti+sj+s = tij and 
ti+5,j = ij,j+5) respectively, while for the latter the same 
relations are the consequence of our phase choices in l|7]). 
With our choice of the orbitals the former matrix ele- 
ments are real while the latter are purely imaginary. The 
matrix elements £55 and £5,10 between the essentially non- 
bonding orbitals are small and are neglected further on. 

Considering the symmetry transformations of the 
structure it is easily found that the 9 distinct non-zero 
first neighbor hoppings are t\g = tj« = t®\ = t 



1 18 
t 10 



-t 



10 



2.10 
10 



- t 01 

- 6 19 

+10 

l 38 



+01 +10 

" c 19j c 27 



10 _ +01 _ +01 
16 — b 16 ~ 
+ 10 _ +01 _ 

l 29 — L 28 ~~ ' ' ' ' 

+01 _ . . . +10 _ +01 _ . . . „ j 

''49 — 5 ^49 — L 38 — ' diiu 

where we did not spell out all re- 
lations which follow rather obviously from the above 
mentioned symmetry transformations of the structure 
and of the orbital products from which also the rela- 
m + n odd follow. As an ex- 



+ 1U +U1 

''4,10 — t 3,10 



tions *™" +5 



j-mn 



+10 

' 18 



-t\s, is ob- 



ample, one of the above relations, t] 
taincd by a 7r-rotation 2 y which transforms (xy) + into 
— (xy) + and (ixz)~ into itself shifted by a translation 
(m,n) = (—2,0). — tfl i s obtained by an 7712- 

reflcction. As an example of the last mentioned type 
(f™ 5 = t%ft_ 5 ), t\l = t\% is obtained by a 7r-rotation 
2 X changing the sign of the orbital (xy) + but not that of 
(— ixz)~, and a subsequent shift by T x mapping (xy) + to 
(xy)~ and (— ixz)~ to {{ixz) + )* . Since (ixz) + couples to 
(xy)~ via V-, the T^-shift changes the sign of the matrix 
element back. On the other hand, in tl° 1Q = t^j no sign 
change of orbitals happens in the 7r-rotation 2 V , and no 
sign change in the shift by T x since this time the orbitals 
couple via the kinetic energy operator and V+. 

There are also 9 distinct non-zero second neighbor 
hoppings tjj,i,j = 1,...,5, while t\\ = 0, t\\ = 0, 

both latter cases by an 7712- 

11 /ll r,y^A ^mn +mn 



+11 - +11 

reflection 



+11 + 11 

6 13' [ 24 



11 



l 23' 



r 33; r 45 



tU and t%" 



j 1 



m 
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even. For instance for ij = 11 the latter relations are ob- 
tained by an mi-reflection, and, for i or j or both equal 
to 3 or 4, by an additional m2-rcncction. As already 
obtained in the general case, t]^ 5 j +5 = tjj. 

Since we try to limit the number of parameters as much 
as possible, and we focus on the low-energy vicinity of the 
Fermi level, we neglect £55 and £5^0 for the Fe-Fe non- 
bonding bands which stay away from the Fermi level by 
more than half an cV. (Note, however, that for the LDA 
or GGA relaxed structure parameters these bands may 
even cross the Fermi level; so far there are nevertheless 
no reasons to assume that this happens in nature.) The 
hybridization parameters tjs, i = 1, ... ,4 however can- 
not be neglected. They strongly influence some Fermi 
velocities. 

The hoppings t 11 happen essentially through the c/p 
site. Due to the large radial extent of the c/p-4p or- 
bitals, nearest neighbor c/p-c/p hopping along the c/p 
tetrahedron edges is strong, and hence the effective Fe- 
Fe hoppings i 20 , i 21 and t 22 are not yet small. In order 
to fit besides the FS radii also the Fermi velocities rea- 
sonably well, they must be taken into account at least 
for the orbitals 1,3,4 and 6,8,9. 

Fig. [9] shows the two-dimensional BZ of the FeSe slab 
of Fig. [SI A notation 



tation is 



ki = k x + k y , fc 2 = -k x + 
is used, and Bloch sums are defined as 

1 

where S; is the site of the basis orbital c 



(9) 



(10) 



^11 



U++ — 
H 12 ~ 

H ++ - 

n l3 — 

= 

H 22 ~ 
H 23 + = 

H ++ - 

H 25 + = 
H 33 



H 34 

H 35 + 
U++ 



ei + 2i}}(cos k\ + cos fc 2 ) + 
+ 2i 2 ?(cos(2fci) + cos(2fc 2 )), 
0, 

2it}3(sinfci — sinfc2), 
2ii}3(sin k\ + sin k 2 ), 
2t]l > {coski — cos £2), 
£2 + 2t\\ (cos ki + cos ki ) , 
2^3 (sin k\ + sin k 2 ), 
2itH(— sin fci + sin k%), 
0, 

e 3 + 2^3 (cos ki + cos fc 2 ) + 
+ 2tf 3 cos(2fc ffi ) + 2t° 3 l cos(2k y ) + 
+ 4t 22 cos^fe^) cos(2fc y ), 

2t34(cOsfci — COSfc^), 

: 2itg5(sin fci + sin k 2 ), 
e 3 + 2^3 (cos fci + cos fc 2 ) + 

:(2k x ) + 2tf 3 cos(2fc. y ) + 
i^fe^cos^fc,,), 
2itgg(sinfci — sinfc2), 
£5, 

ij 



(11) 



■2t° 2 cos 



k 2 




Figure 9: BZ of the structure of Fig.[8j dashed for the 
Fe structure alone (one Fe atom per cell) and full lines 

for the full FeSe structure. The wave number 
components are scaled so that the symmetry points are 
(k x ,k y ) = (0,0) for T, (7r/2,7r/2) for X, (tt,0) for M and 
(tt, tt) for T', as well as (fci, fc 2 ) = (0, 0) for T, (vr, 0) for 
X, (7r,7r) for M and (2tt,0) for V. 

The Hamiltonian matrix H ++ in Bloch state represen- 



Of course, H ++ is Hcrmitian, and in agreement with 
Hermiticity and our phase choice in one finds that 
tjj is real for tjj = t^j and imaginary for tjj = — tj*. 

Hence, one finds that is real. Together with (|5|) 

this also implies 



H— = H ++ . (12) 



This reality is the consequence of absence of hoppings in 
i?3 direction and does not hold in the more general case. 



8 



For H + one finds 

H^q~ = 2tJg(cos k x + cos k y ) + 

+ 2t\ 6 ( (cos fci + cos fc 2 ) (cos fca; + cos k y ) 

— sin k\ (sin fc^ + sin k y ) + 
+ sin fc 2 (sin /c x — sin fc y ) ) , 

H+~ = 0, 

H^g = 2zi^g sin fc^ , 

Jf-^g — 2z£-^g sin . 

#1,10 = °) 

#2^7 _ = 2^7 (COS fc^ + COS 
7^28 = — 2z't 2 9 sni ^y; 

= 2z^29sinfc ;c , 

-^-2.10 ~ ^2,10 ^ x COS fey), 

.f/^g — — 2i^g cos ~f~ 2i^g cos fcy — I - 

+ 2tgg((cos k\ + cos /C2) cos k x — 

— (sin k\ — sin fc 2 ) sin k x ) + 
+ 2^4g((cos k\ + cos £2) cos k y — 

— (sin &i + sin ki ) sin k y ) , 

H 3 +- = 0, 

^3,10 = 2zi4 10 sin 

ff^Q — 2i^g cos /c^ H- ^^38 cos /cy H- 

+ 2£4g((cos k\ + cos /C2) cos k x — 

— (sinfci — sinfc2)sinfc K ) + 
+ 2fgg((cos k\ + cos /C2) cos k y — 

— (sin fci + sin ) sin /c y ) , 
-ff^io = 2it^ iq sinfca;, 

<fo = 0. 

is again real, and hence 

H+~=H-+. 
In total, the block structure (0 specializes to 




(13) 



(14) 



(15) 



This Hamiltonian is easily block diagonalized. Con- 
sider the Bloch sums 



(16) 



i = 1, . . . ,5, 



with 4>f from ^ and form the combinations 



where s stands for symmetric and a for alternating. Over- 
lap of neighboring local orbitals is usually neglected in a 
tb approach. 

If one, more generally, denotes with Wannier func- 
tions with the same crystal site symmetry as then 
overlap in the original basis ([?]) is automatically included. 
This is done in all what follows. (See also next section.) 

Then, fl7|) immediately yields 



H B8 = H ++ + H +- 

H sa = = H as . 



jjaa = H ++ _ jj-i 



(18) 



The next section deals with the determination of tb 
parameters from FPL08 band structure results. Here, 
the tb model jHHIBl) for the case of FeSe is analyzed as 
a typical example. The corresponding tb parameters are 







t 11 - 
'11 — 


0.086, 


,10 

' 16 


= -0.063, 






,20 

'ii — 


-0.028, 


' lfi 


= 0.017, 






t 11 - 

'13 — 


-0.056*, 


.10 
'18 


= 0.305i, 






t 11 - 
'15 — 


-0.109, 


'27 


= -0.412, 


ei = 


0.014, 


t 11 - 
'22 — 


-0.066, 


,10 
'29 


= -0.364i, 


£2 = 


-0.539, 


t 11 - 
'23 ~ 


0.089i, 


,10 

I 2,10 


= 0.338, 


£3 = 


0.020, 


t 11 - 
'33 — 


0.232, 


,10 

'38 


= 0.080, 


e.5 = 


-0.581 


.20 

'33 — 


0.009, 


,21 

'38 


= 0.016, 






.02 

'33 — 


-0.045 


.10 
'49 


= 0.311, 






,22 
'33 


0.027, 


c 49 


= -0.019, 






< _ 

'34 — 


0.099, 


.10 

'4,10 


= 0.180i, 






t 11 - 
'35 — 


0.146* 




all in eV. 



(19) 



h+i -I- A-^ 



(17) 



Only in a few cases the on-site energies e, and the hop- 
ping matrix elements t directly obtained from the Wan- 
nier approach described in the next section are slightly 
renormalized in order to grossly account for the neglected 
hoppings to higher neighbors. 

Fig. [10] shows the obtained 2D two to five neighbor tb 
bands in comparison with the full LDA bands of the 3D 
FeSe crystal. For good reasons the focus of this treatment 
is on the vicinity of the Fermi level. On the first glance 
it seems not to be an optimal fit there, in particular not 
around the points M and A (see Fig. [5] for M, the points 
Z, R, A are at k z = ir above T, X, M). However, due 
to hybridization with the Se-4p states the Fe-Fe anti- 
bonding bands have a sizable fc 2 -dispcrsion which cannot 
be provided by a 2D tb model. The quality on average 
of neglect of fc 2 -dispersion of this tb model is also seen 
from the comparison of the DOS in Fig. fTTI showing that 
the fit is not unreasonable close to the Fermi level. 

A further feature of the 2D tb model is the block diag- 
onalization (|18p of the Hamiltonian matrix by means of 
lattice symmetry. There are two five-band groups in this 
model which do not interact in the whole BZ. One group 
comes from symmetric binary Wannier states (|17j) and 
one group from alternating ones. From what was said in 
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too, for which there is no continuous outfolding. There 
is no additional benefit besides that. 



Figure 10: (color online) Comparison of the 2D tb 
bands with the full 3D LDA bands of FeSe. 





-1 

energy in eV 

Figure 11: (color online) Comparison of the 2D tb DOS 
with the full 3D LDA DOS of FeSe. 



Figure 12: (color online) The ss- and aa-band groups 
corresponding to the Hamiltonian blocks (p~6|) for FeSe. 

In summary so far. without considerably worsening the 
fit. no less than the 27 tb parameters (fl"9)) will do in 
reasonably modeling the electronic structure of the iron 
chalcogenides/pnictides, even if one focusses on the Fermi 
level only. This is, unfortunately in stark contrast to the 
situation for the cuprates. 

As can be inferred from Fig. [TU1 a 2D tb fit for 
FeSe cannot be very satisfactory for quantitative stud- 
ies even in an energy window of ±0.1 eV around the 
Fermi level, relevant for low temperature many-body ap- 
proaches. This is different for the families 1111 and 122. 
In this respect, although their crystal structures are more 
complex their low-energy electronic structure is simpler 
and really distinctly 2D. 

For LaOFeAs, the 2D tb parameters are obtained as 



the Introduction it follows that each of the two groups 
contains two Fe-Fe bonding and two Fe-Fe anti-bonding 
bands as well as one essentially non-bonding band. The 
two groups are shown on Fig. [T2l Within each group 
there is no further decoupling in the whole BZ. In other 
words, a further reduction of the number of bands of a 
tb model necessarily demands inclusion of (many) higher 
than second neighbors to mimic the fc-dependent inter- 
action with the omitted bands. 

On the side faces of the BZ the two groups are degener- 
ate. As was demonstrated in Section II, this degeneracy 
is a general property of the P4/nmm space group. It en- 
ables to fold out smoothly the ten bands into five bands 
in a doubled BZ (Jones-zone^). One may chose the ss- 
bands in the BZ proper and the aa-bands folded out, or 
vice versa. It is now clearly seen that in any case this is 
not an out-folding of anti-bonding bands against bond- 
ing bands. Were it so, a pseudo-gap between bonding 
and anti-bonding bands could not develop. 

The other simplification, the reduction of a 10 x 10 
Hamiltonian matrix to a matrix structure with two 5x5 
blocks H ++ (Hcrmitian) and H A (symmetric), is al- 
ready guaranteed by (4) which holds for the 122 family 



ei = 
£2 = 
e.3 = 
e.5 = 



0.163, 
-0.407, 

0.053, 
-0.196 



ii 

20 
11 

11 
13 
11 
15 
11 
22 

11 

23 

11 

20 
:i:S 

02 
:i:S 

22 
33 
11 
34 



l 35 — 



0.120, 
-0.029, 
-0.01.4i, 
-0.172, 
-0.038, 
0.079i, 
0.235, 
0.023, 
-0.025, 
0.032, 
0.094, 
O.llli 



,1U _ 

r 16 — 



,10 _ 
t 2,10 ~~ 

,10 _ 



1 38 
,21 

c 38 

,10 
c 49 

,21 
c 49 



,10 _ 

c 4,10 — 



-0.167, 
0.027, 
0.224i, 

-0.348, 

-0.315i, 
0.296, 
0.093, 
0.026, 
0.335, 

-0.008, 
0.126i, 



all in eV. 



(20) 



The corresponding low-energy 2D bands are compared 
with the full 3D LDA bands in Fig. [T31 Since here the 
FeAs layers are separated by LaO layers which have no 
states in the vicinity of the Fermi level, the full band 
structure has a pronounced 2D character there, and con- 
sequently the 2D tb approximation is of high quality. It 



10 



does not need improvement by considering 3D disper- 
sion. Nevertheless, for the LDA or GGA relaxed struc- 
ture instead of the here considered experimental Wyckoff 
positions one of the z 2 -derived bands crosses the Fermi 
level instead of the xy-derived one. This is why many 
published band structures show a hole pocket around Z 
instead of the third hole cylinder. 



On the expense of an increase of the number of tb 
parameters to 44, of course a better overall fit of the 
2D band structure in a larger energy window can be 
obtained^ 



The 2D tb parameters for BaFe2As2 are obtained as 




Figure 13: (color online) Comparison of the 2D tb 
bands with the full 3D LDA bands for LaOFcAs. (top) 
full energy window, (bottom) zoom in the low energy 
region. 



Somewhat surprisingly, the situation for BaFe2As2 in 
a narrower energy window is similar although not equally 
perfect. This allows to treat the low-energy band struc- 
ture of the 122 family within the same space group 
P4/nmm (and BZ) which is the space group for one triple 
layer of FeAs, although the space group of the full 3D 
crystal is I4/mmm (due to an ro-refiection of adjacently 
stacked FeAs triple layers) . As a convenient consequence, 
the Hamiltonian (fTTl [T3"l) with its simplifying block struc- 
ture (|15|) applies. 







t 11 - 


0.135, 


.10 

' us 


= -0.196, 






,20 
l ll ~ 


-0.027, 


'16 


= 0.042, 






t 11 - 


-0.024z, 


.10 
' 18 


= 0.218i, 






t 11 - 
l 15 — 


-0.131, 


.10 
'27 


= -0.355, 


ei = 


ii. 1 72. 


t U - 
L 22 — 


-0.131, 


,10 
l 2<) 


= -0.365*, 


(2 = 


-0.236, 


t 11 - 
c 23 — 


0.103i, 


.10 

I 2,10 


= 0.265, 


e.3 = 


0.000, 


l 33 — 


0.204, 


,10 

c 38 


= 0.065, 


e.5 


-0.590 


.20 

l 33 ~ 


0.034, 


/ 21 

L 38 


= 0.020, 






.02 

''33 — 


-0.048, 


.10 
'40 


= 0.312, 






.22 

c 33 — 


0.024, 


t 21 
'4!) 


= -0.024, 






f 11 - 
\34 — 


0.118, 


.10 

I 4,10 


= 0.080i, 






t 11 - 
E 35 — 


0.078i 




all in eV, 



(21) 



and the low-energy bands together with the full 3D LDA 
bands are shown on Fig. [TU (Note, that we used the 
same path along the BZ as for the other cases, although 
the symmetry of the 122 system differs. This, however, 
facilitates comparision.) Of course, in both Figs. [131 and 




Figure 14: (color online) Like Fig. [12] for BaFe2As2. 
Above the Fermi level the full LDA bands are strongly 
hybridized with Ba 5d states. 



[14]the tb bands on lines T-X-M-r and Z-R-A-Z are iden- 
tical, but the full 3D LDA bands arc not and the plot in 
particular also shows the distinct 2D character of the full 
bands and the quality of the 2D fit. 

For later use we also give the 2D tb parameters for 
LiFcAs valid in connection with the Hamiltonian (fTTl [T3"l) 
although without a 3D generalization provided in Sec. V 
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the fit would not even qualitatively be correct. They are 







+ 11 

Hi = 


U.U ( y, 


,10 
' 16 


n ai £ 
= — U.UlD, 






hi = 


a aoa 
0.020, 


,21 
' 16 


= 0.013, 






,n 

'13 — 


-0.090i, 


, 1 
r 

' 18 


= 0.281?, 






*1B = 


—0.060, 


,10 
'27 


= -0.404, 


ei = 


-0.188, 


,11 

'22 — 


a no i 
— O.Ooz, 


,10 
'20 


= — 0.353l, 


£2 = 


—0.521, 


+ 11 

t 23 — 


U.Uo it, 


,10 

''2,10 


— U.olo, 


£3 = 


U.ZUU, 


,11 

'33 ~ 


A i7r; 
U.zVo, 


,10 

'38 


= U.lzo, 


e 5 = 


— n fifiQ 

u.uuy 


,20 

'33 — 


-0.002, 


,21 

'38 


= 0.056, 






,02 

'33 — 


-0.107, 


,10 
'40 


= 0.359, 






,22 

'33 — 


0.012, 


t 21 
'4!) 


= -0.048, 






'34 — 


0.102, 


,10 

I 4,10 


= 0.190i, 






t 11 - 
'35 — 


0.136« 




all in eV. 
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IV. FPLO8 WANNIER FUNCTIONS 

The determination of a Wannier function (WF) basis 
of states \Ri) with function index i at center Si in unit 
cell R, though formally being a simple Fourier transform 
of the Bloch states \kv), 

\m) = Y / \^)u!; l c-^ R , (23) 

h.u 

becomes an art due to the inherent ambiguity of the 
choice of the band index and fc-vector dependent phase 
factors U^. 

There are basically three goals pursued, when calcu- 
lating WFs, the first being to obtain functions, which 
form a suitable basis for the setup of model Hamiltonians. 
This aim usually implies that the resulting WFs should 
resemble basis functions, which have a certain chemical 
characteristic. The second goal is to reduce the degrees 
of freedom in the resulting model by projecting out high 
energy sectors resulting in few-band models. The last 
goal is to obtain model Hamiltonians with as few as pos- 
sible parameters, which is essentially the requirement of 
maximal localization of the WFs. RefJ^ introduced a 
general method of calculating maximally localized WFs. 
Due to the nature of the problem at hand this approach 
is rather involved. 

In order to fulfill the condition on the chemical char- 
acteristics it is in general desirable to pose symmetry 
restrictions on the WFs. In a local orbital scheme with 
an optimized chemical basis set the basis orbitals usually 
show those characteristics, however, the resulting basis 
is non-orthogonal. Nevertheless, the local basis orbitals 
can be used as a starting point for the definition of WFs 
of a certain symmetry, which by construction have a high 
degree of localization. 

We follow the ideas sketched in Ref^ and introduce 
a projection of the Bloch functions \kv) onto local test 



functions \Si) . The resulting matrix elements (kv\Si) 
enter the phase factors U^. They assign weights to the 
Bloch states entering the Fourier transform Eq. (|2"3"]) such 
that the bands with the largest test function character 
contribute the most. It is very conceivable that opti- 
mized local basis orbitals are among the best possible 
choices of test functions. Thus in the simplest case we 
choose \Si) to be a local orbital at the site S with orbital 
quantum numbers i. Alternatively, a linear combination 
of local orbitals (molecular orbital, MO) centered at the 
WF center S could be used, in which case i denotes the 
characteristics of the MO. 

The choice of the projector singles out all the bands ex- 
hibiting a particular character. In order to construct few- 
band models it is often desirable to additionally project 
out certain parts of the band complex with a particular 
character, say one is interested in WFs describing the 
anti-bonding part only. This is achieved by introducing 
energy windows described by functions hgi{tku), which 
are basically unity in the relevant energy range and zero 
everywhere else. 

All the projectors put together define the Bloch sums 
corresponding to the WFs 

\ki} = \ku)h S i{e h u){kv\Si) (24) 

which in general will be non-orthogonal. A subsequent 
symmetric orthogonalization defines the orthogonal WF 
Bloch sums. 

\M) = j2m(s- l/ % (25) 

j 

where Sjfc = {kj\ki}. The resulting phase factors 

U5t = E M^) {kv\Sj) (S- 1 / 2 )* (26) 
3 

are then used in the Fourier transform Eq. (f2"3"|) to calcu- 
late the WFs. They have predominately the character of 
the test functions | Si) and approximate the bands which 
exhibit this character and lay in the chosen energy win- 
dow. The WFs will transform according to the symmetry 
of the test functions. Due to the construction the largest 
part of the WF is actually formed by the test function 
itself. Since, we use optimized chemical orbitals as test 
functions, the resulting WF will have a high degree of 
localization. 

Furthermore, this construction results in WF band 
structures which do not only fit the band energies but 
also the orbital character of the bands. A band, which 
has predominantly a particular orbital character in the 
LDA bands will have the character of the corresponding 
Wannier function in the WF bands. Since the WFs are 
constructed to resemble the orbitals the Wannier charac- 
ter reflects the original orbital character. We confirmed 
this behaviour explicitely by comparing orbital/WF pro- 
jected bands. 

In the present work, we used the Fe 3d orbitals as test 
functions. The energy window was chosen such that the 
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c/p-4p bonding bands where excluded from the energy 
window and only the 10 Fe-dominated bands were in- 
cluded. 



V. THREE-DIMENSIONAL DISPERSION 

As indicated in Sec. Ill, there are reasons to consider 
the /c z -dispersion close to the Fermi level in the family 
11, and this becomes mandatory for the family 111. This 
increases the number of parameters considerably since it 
is connected with hopping through several layers, and, 
since there are no vertical bonds, each of these hoppings 
adds also one vector ±Ui/2 or ±i?2/2 (for hops from Fe 



to c/p) or ±T X or ±T y (for hops from c/p to c/p or to a 
cation) . 

In the 11 and 111 families, at least three hops are in- 
volved from an Fe layer to the next Fe layer (t rst with 
t = ±1), and hence t rs±1 with (r,s) = (0,0) or (±1,0) 
or (0,±1) or (±1,±1) or (±2,0) or (0,±2) or (±1,±2) 
or (±2, ±1) are of the same order of magnitude and have 
to be considered together. (Li or Na orbitals of the 111 
family are not involved in the considered case.) For the 
vicinity of the Fermi level it suffices to consider the k z - 
dispcrsion of the xy, xz and yz bands. The additions to 
the Hamiltonian matrix (jlll I13[) are (k z = tt at BZ-point 
Z) 



H ll 


— W++ 


+ (2i°? 1 + 4^} 1 (cosfc 1 4 


cos fe) 


+ itfl 1 (cos k x + cos k y )) 


W+ + 
H 13 


U+ + 
~ H 13 


-4t\f sin(2fc. y )sin k z , 








- H 1A 


- itfl 1 sin(2k x ) sink z , 








- H ++ 

— -"33 


+ fit™ 1 + At™ 1 cos(2k x ) 




cos(2k y )) cosfc z , 




- H ++ 


+ (24f + AtH 1 cos(2/s x ) 




cos(2fc y )) cos k z , 



H 19 
#3 + 8~ 



#39 



H 



49 



H 16 



-4^ (cos k x - 
-2^i(( cos ( fcl 



H+T =H+--4(tl° 8 1 sink 



- cos k y ) cos k z + 

- cos(fci + k x )) exp(ik z ) + (cos(fc : 

^ „• I. , o-j-121 l-n. , \ 



#19 

# 3 V 



#39 



+ ky)- 

- 1^ 1 sinfcj,) sinfcz 

- 4(t^g 1 sin k x + t^ 1 sin k y ) sin k z - 

- 4(tgg 1 cos k x + t]^ 1 cos k y ) cos k z 

- iitlg 1 (cos k x + cos k y ) sin k z , 



-2zt}g 1 (sin(fci 



'2 T rv y 

ky) exp(i/v z 
ik 



2it}g 1 (sin(fci + k x ) exp 

k x )exp(ik 



- 2t^ 1 (cos(fc 1 

" 2t 4 g 



(cos(fci + k y ) exp(ik z 

4(t4g 1 cos k x + tl^ 1 cos k y ) cos k z + 2t 4 g (cos(fci + k x ) cxp(ik z 

(cos(fci + k y ) cxp(ifc 2 



k y ) + cos(fc 2 - k x ) exp(-ik z )), 
k z ) - sin(fc 2 + ky) * exp(— ik z )), 
sin(fc 2 - k x ) * cxp(-ifc z )) , 

- cos(fc2 — k x ) cxp(—ik z )) + 

- cos(fc 2 + ky) exp(—ik z )) , 



9.121 

4- 2t 121 



- COs(fc2 

- COs(fc2 



- k x ) exp(-ik z )) ^ 

- k y ) exp(-ik z )) . 



(27) 



Since the Hamiltonian matrix is now complex, (fTS"]) does cV) 
not hold any more, but (0 [6]) is still true. 



C= 0, 

,oo! = n ^ = -0.017, 

" " ' t\ 01 = 0.009*, 

t 111 =0 

tf 11 = 0.017, " 

II, 41 1 = 0.03H, 

t^= 0.030*, |° (28) 

1 tl 01 = 0.006, 

4^= 0.011, ™ 

om 41 = -0.003, 

tH 1 = -0.008, ™, 

tlf = 0.020, f m 

t™ 1 = 0.025, 

tH 1 = 0- 006 - 



For FeSc, of the 17 new tb parameters contained in The comparison of the corresponding low-energy tb 
(|2"T|) the 14 non-negligible parameters are (again all in bands with the full LDA bands is shown on Fig. [T5J The 
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fit close to the Fermi level is probably about the best 
which can be achieved in this case with a total of 41 pa- 
rameters. The 10 Fe-related bands of FeSe are separated 
from the rest of the bandstructure by gaps (see Figs. [5] 
and [6]) . If one takes the full set of several hundred of tb 
parameters calculated from our WFs, the red lines merge 
nearly perfectly the black ones in Fig. [TO] for all 10 Fe-3d 
related bands shown. This has been checked. 
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Figure 15: (color online) Comparison of the 
fc 2 -dispersed tb bands with the full LDA bands for 
FeSe. (top) full energy window, (bottom) zoom in the 
low energy region. 



bands in Fig. [TOJ 

Besides the Fermi radii around M and A now being 
markedly different, there now appears a Fermi pocked 
around T with a FS crossing even the line L-Z. One radius 
around T could be improved on the expense of many more 
hopping parameters over larger distances only, which is 
again due to the far extension of the As-Ap orbitals. 




Figure 16: (color online) Like Fig. [15] for LiFeAs. 



For LiFeAs, 16 of the 17 new parameters contained in 
(|2"T)) are relevant. They are 



ooi _ 

n — 

in _ 

n — 

201 _ 
11 — 

201 _ 

14 — 

001 _ 
33 — 

+201 _ 
L 33 — 
+021 _ 
L 33 



0.070. 
0.020, 
0.005, 
0.025i, 
-0.004, 
-0.003, 
0.031, 



101 

16 

211 

16 

101 

18 

101 
19 
211 
19 

101 

38 

211 

38 

101 

39 

101 

49 

211 

49 



-0.018, 
-0.025, 
0.008i, 
0.020i, 
0.022i, 
0, 

-0.014, 
0.017, 
0.016, 
0.043. 



(29) 



Let us finally shortly re-examine the 122 system. With 
a few z-hopping matrix elements, the tb fit in the wide- 
energy window shown in Fig. 14 cannot really be im- 
proved since the Ba 5d states, strongly hybridizing right 
above Fermi level, are very extented. Nevertheless, a few 
of those hoppings yield a better overall picture. The ma- 
trix elements in addition to (21) are 



,101 

l n 

+211 
L 14 

+ 101 

z 33 



0.061, i^ 1 = 
-0.029i, t™ 1 = 



0.030, 



+m - 
He — 



-0.161, 
0.080, 
-0.041. 



(30) 



The corresponding bands are compared with the full LDA 



Due to the different structure the addings to the Hamil- 
tonian matrix are different from (27). They are 
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= H{{ + (cos k x + cos k y ) cos k z , 
Hi~3~ = + 2t^4 1 (i sin k x cos(2fc y ) cos k z + cos k x sm{2k y ) sin fc z ), 

= H^±~ + 2ti\ 1 (sin(2k x ) cos k y sin fc z + i cos(2fc.,;) sin k y cos fc z ), 

= #|| + 443 1 cos k x cos fc z , 
#/ 4 + = Hll + it™ 1 cos k y cos fc z , 

-^icT = -^ieT + ^t^g" 1 cos fc z + 2tjg (cos fci exp(— ifc z ) + cos fc 2 exp(ik z )), 

-^38 = -^38 + ^t^g 1 COS fc z , 

#3 + 8 " =H+T + 2t° 8 ll cosfc z . 



This 3D tb fit is shown on Fig. \T7\ As is seen, there is 
one band (essentially of Ba 5d character) comming down 
close to the Fermi level which is not represented by the 
Fe-3d tb fit. At variance to the case of LiFeAs, here the 
3D dispersing bands are, however, all unoccupied (and 
even move away from Fermi level under hole doping). 




Figure 17: (color online) Comparison of a k z dispersed 
33 parameter tb fit of the Fe 3d bands with the full 
LDA bands for BaFe 2 As 2 . 

If one plays with the c/p-Wyckoff parameter or with 
larger ordered magnetic moments causing large addi- 
tional exchange potentials and band shifts, the 3D be- 
havior and maybe a FS crossing the line T-Z can reap- 
pear even for the 1111 and 122 families. This would put 
serious questions on the practicability of the tb approach 
at all, and one would probably have to go back to a full 
numerical WF treatment of the band structured 



VI. CONCLUSIONS 

We have demonstrated that the band structures of the 
non-magnetic undoped iron based superconductors for all 
four families have a quite complex multi-orbital character 
with all 10 Fe-3d orbitals per unit cell involved in band 
states in a 0.1 eV vicinity of the Fermi level. Thereby, 
the states at the Fermi level are mainly of xz, yz and 
xy character, respectively. However, departing from the 
Fermi level, the other Fe-3ci orbitals and even the c/p-4p 
orbitals start to hybridize and largely influence the Fermi 
velocities. 

While the most investigated families 1111 and 122 
exhibit clear 2D behavior in the low electronic energy 
regime, this is much less the case for the family 11 and 
not at all true for the family 111. Minimum quantita- 
tively correct tb models arc presented for FcSc, LiFeAs, 
LaOFeAs and BaFe2As2 as representatives for the four 
families. 

A minimum number of 27 tb parameters seems not 
to be reducible without severe loss of accuracy of the 
corresponding Fermi radii, Fermi velocities and orbital 
character on the FS. For the families 11 and 111 even 
that is not sufficient, and the minimum number of tb 
parameters goes beyond 40. 
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